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1. Introduction 

Random walks are ubiquitous in fields like physics, chemistry, biology, and economics. 
In a generic version of a random walk with discrete steps, every time a walker arrives 
at a crossroad, he has to choose the route to take. After several crossings and choices, 
made for example by flipping a coin [with heads (tails) leading to a step to the left 
(right)], he will have followed one out of many possible paths. For a quantum walker, 
in contrast, the result of each coin toss is a superposition of heads and tails. That is, 
the quantum coin takes both states simultaneously and therefore the walker follows all 
the possible paths. pQ EJ El H [5] 

The results of a classical walk are intuitive and relatively easy to comprehend. 
Quantum random walks, however, show strange unintuitive results due to the possibility 
for the walker to meet himself, and either increase his probability to be at this crossing 
or even disappear. There is an increasing interest to understand quantum counterparts 
of classical random walks and to find out whether quantum effects are already exploited 
by nature. It has been suggested that single photon random walks have applications 
in quantum computing[6], [7]. Additionally, quantum random walks with two identical 
particles have been shown to solve the graph isomorphism problem[8]. However, previous 
experimental realizations have been limited to single-particle quantum walks, which have 
an exact mapping to classical wave phenomena [3], and therefore cannot provide any 
advantage (computational or otherwise) as a result of uniquely quantum mechanical 
behavior. 

Quantum random walks with photons can be achieved by using a staggered array 
of 50/50 beam splitters each representing a step left or right. Now two walkers realized 
as photons were observed to pass through the same optical path network. For the first 
time, the two indistinguishable walkers were shown to interfere with each other [I]. In 
contrast, for quantum walks of more than one indistinguishable particle, classical theory 
no longer provides a sufficient description. Quantum theory predicts that probability 
amplitudes interfere, leading to distinctly non-classical correlations [Hj. 

In this paper, we show that there are interesting aspects to these higher level multi- 
photon random walks. Dirac's dictum, that every photon only interferes with itself is 
wrong. Hong, Ou, and Mandel showed that two photons incident on a beam splitter 
tend to interfere and both go one way and both go the other way [9]. For higher photon 
numbers, a generalized Hong- Ou- Mandel effect prevails. However, the outcomes of these 
higher level multi-photon random walks are harder to visualize. We investigate onefold, 
twofold, and threefold probability coincidence detections in order to extract quantum 
information about the walks with multiple walkers from the output of these random 
walks. Based on "Feyman diagrams", we show that although quantum random walks 
run in parallel, coincidence detection at a particular site responds to specific paths only. 

In future research, we plan to explore higher order correlations, use more generalized 
number resolving detection schemes, and investigate genetic algorithms in order to steer 
photons using phase shifters at key points within the array to create desired outputs. 
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2. Quantum random walks 



2.1. Model 



A peg board and a ping-pong ball, as shown in Figure 1(a), is an example of a 



classical random walk, where the ping-pong ball is the walker that bounces left or 
right off the pegs on the board. Repeated walks through this setup result in outcomes 
that are distributed according to the Gaussian probability distribution with maximum 



probability right beneath the input peg as shown in Figure 1(b) 




Position 



(a) 



(b) 



Figure 1. (a) Classical random walk using a peg board and tennis balls, (b) Result 
of a classical random walk is a Gaussian distribution. The maximum of the Gaussian 
corresponds to the location directly below the starting position. 

In the case of quantum random walks [H [2], we replace the pegs with 50/50 beam 
splitters and the ping-pong balls with photons. From a classical stand point they 
operate by splitting incident light so that half of the incident energy is reflected with 
an additional 7r/2 phase shift acquired by the field from the beam splitter, and half is 



transmitted through the beam splitter (see Figure 2(a)). 



There are many ways to model a beam splitter quantum mechanically. Yet, we 
have chosen to use the symmetrical model, which follows the physical properties of 
laboratory, optical, beam splitters, instead of an asymmetric, non-physical Hadamard 
model. This way, our results can be easily reproduced in an all-optical experiment. A 
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Figure 2. (a) Classical description of a beam splitter; (b) Quantum mechanical 
description of a beam splitter with the input and output field operators shown. 
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Figure 3. Pyramid structure of quantum random walk model using beam-splitters. 



quantum-mechanical description, adopted here, of such an experimental beam splitter 



is given by the transformation for the field operators illustrated in Figure 2(b) , [TU] 



a in = 4fKut + ^out)' 

f (1) 

6 in = 7f( a out+ i6 out)- 

The terms that contain the factor of i come from the fact that the incident photon is 
reflected off the beam-splitter instead of being transmitted through it, which gives the 
photon an additional 7r/2 relative phase [TU]. Here individual photons after interactions 
with a beam-splitter end up in a superposition of being simultaneously transmitted and 
reflected. Consequently, the probability distributions for the outcomes are no longer 
Gaussian, but they multiphoton demonstrate quantum interference effects instead. 

We start with a two-mode quantum state of light with N photons at the left port 
of a beam-splitter and M at the right, written, \N,M) = ^^ lM , a ^ N b' M 1 0, 0) . Here the 
input photons correspond to the creation operators, which after the first beam-splitter 
become 



1 / ia) + \ 1 / at + i&t 



/ iVT V a/2 7 y/W\ \ V2 



(2) 



where creation operators a) and W corresponds to the photon exiting the left and 
right output port, respectively, of the beam-splitter [10]. This process is cascaded 
down through additional beam-splitters into a pyramid shape as shown in Figure |3j 
Hence, in order to propagate photons through the system, we apply the beam-splitter 
transformation (TlJ recursively onto a) and b' , while adding new creation operators that 
correspond to photons exiting from output ports of these beam-splitters, until we have 
reached the desired level. At this point, we turn the result from operator into the state 
notation by acting with the creation operator on the vacuum state, a) n |0) = y/n\\n). 



2.2. Detection of Outcomes: Correlation Functions 

The outcome of classical random walk is a walker exiting at a certain port. In a quantum 
random walk, the outcome is a state that describes the statistics of possible outcomes. In 
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order to recover this state, one can use conventional photodiodes. However, they do not 
distinguish the number of photons that actually hit them. In particular, the signal from 
a detector at position m is proportional to the mean photon number (^l a|„a m \ip). In 
single-photon experiments, the mean photon number is equal to the onefold probability 
that detector m receives a photon. The probability distribution over all detectors 
provides all the information. It is important to note that single-photon and multi-photon 
walks are described by identical probability distributions. In experiments with multiple 
quantum walkers however, the onefold detection scheme throws away information, as it 
does not measure the exact number of photons that impact the detector and ignores 
what happens to other detectors. 

Quantum random walks with multiple photons have unambiguous differences in 
the characteristics of the outcomes due to the interference between the two or more 
indistinguishable walkers. This new property of two, indistinguishable, interfering 
photons cannot be revealed with onefold probability. The twofold probability expands 
upon the onefold probability, where the twofold probability at detectors m and n is equal 
to (^| d^d^dndm \ip) [121 E321 H] • The twofold probability reveals the correlations between 
simultaneous events at detectors m and n or acknowledges that two photons were 
received by the same detector. The twofold probability detection scheme is described 
by the joint probability that two detectors fire at the same time. This means that the 
twofold joint probability at detectors m and n represents the probability that detectors 
m and n each measure at least one photon simultaneously. Due to symmetry, the twofold 
probability at detectors m and n is equal to the twofold probability at detectors n and 
m. 

In the case where we have detector m = n, the twofold probability at detectors 
m and m represents the joint probability that detector m fires and detector m fires. 
We interpret this as detector m receiving exactly two photons. After dividing by two 
due to the bosonic nature of photons, this interpretation holds well for two walkers and 
the twofold probability can be measured with conventional photo-diodes by repeating 
the same walk multiple numbers of times and counting the frequencies of n, m events. 
Similar to the onefold probability, twofold probability throws away any information that 
involves three or more detectors firing simultaneously and it also can not distinguish 
between the cases where two photons impact a detector, or the cases where greater than 
two photons impact a detector. 

Onefold and twofold probabilities have been thoroughly studied previously, [H 
[TTl [T2l [13] and are mentioned in this paper for completeness. It is only logical 
for multi-walker quantum random walks to consider higher-order correlations. We 
present a case of threefold probability, which has not been previously studied, and 
discuss the visualization of the three dimensional results. It is defined similarly to 
lower-order probabilities; threefold probability at detectors m, n, and / is equal to 
{^\d\ n d\jd\did n d m \'ijj) . The threefold probability at detectors m, n, and I represents 
the probability of measuring a photon at detectors m, n, and / simultaneously. The 
threefold probability also has symmetry similar to the twofold probability, as the 
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threefold probability at detectors m, n, and / is equal to the threefold probability at 
detectors n, I, and m, and these are equal to any other permutations of m, n, and I as 
well. 

3. Results 

3.1. Onefold probability 

As it was pointed out in Section [2^21 onefold probabilities with more than one photon in 
the system produce similar results without providing new information. Because of this, 
we only look at input states with one photon, |1, 0). A flipped input state, |0, 1), results 
in a flipped output state, due to our choice for symmetrical beam-splitters. Figure 4(a) 
shows the onefold probabilities of a |1, 0) input state propagated through three levels of 
beam-splitters. This probability distribution is different from the Gaussian distribution 
obtained in classical walks. We can see that even with a single photon, constructive and 
destructive interference still occurs. This interference is caused by the photon taking 
multiple paths to arrive at a specific detector. 

3.2. Twofold probability 

Quantum random walks with multiple walkers require higher order correlation functions 
to reveal additional interference effects that take place. For a two walker random walk, 



|2, 0) input state produces results similar to |1, 0). Figure 5(a) shows a plot of our twofold 



probabilities for a |1, 1) input state after three levels of propagation. The horizontal and 
vertical axes of the plot correspond to the two detectors m and n. The symmetry can 
be seen down the diagonal of the plot. 

3.3. Threefold probability 

Increasing the number of walkers leads to higher order probability functions being 
deployed to reveal higher order interference effects that are a result of higher number 
of paths traveled by multiple walkers instead of individual paths traveled by a single 
walker. If we consider the case of three walkers, to visualize the results of the threefold 
probabilities on a specific input state, we must find some way of displaying probabilities 
at each of the three detectors independently. We overcome with this challenge by 



essentially stacking plots similar to our twofold probability plots in Figure 5(a), each 
one keeping the third detector constant, and leaving space between each one in order 
to see into it. We call this the "exploded Rubik's cube" plot. Each of the three spatial 
dimensions are used to represent the three detectors m, n, and I, and the color of each 
square represents the threefold probability at those detectors. An example of this type 



of plot is shown in Figure 6(c). Figure 6(a) and Figure 6(b) show the symmetry along 



the space diagonal that is present in all threefold probability plots. Using these plots 
we can also easily see interference effects that carry over from the onefold and twofold 
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Figure 4. (a) Plot of onefold probabilities showing single photon interference of |1, 0) 
input state propagated through three levels of beam-splitters. Complete destructive 
interference is shown at detector D3 while constructive interference occurs at detector 
D4. Two Fcynman diagrams required to explain the complete destructive interference 
at detector D3 (b,c). Feynman diagrams explaining the constructive interference at 
detector D4 (d,e). Only one path is shown per diagram since there is only one photon 
in the system. 



probabilities. For example, if for a specific input state and propagation level, the onefold 
probability at detector a is zero, then for the same input state, the threefold probabilities 
at detectors m, n, and I will equal zero for all m, n, or I = a. You can easily see this 
effect in the threefold probability plots, as entire planes will be zero. Similarly, if the 
twofold probability for a combination of detectors is, then an entire row and column in 
the threefold probabilities will be zero for that same combination of detectors. 



4. Discussion 



4-1. Feynman Diagrams 



Coincidence detections with conventional photodiodes requires statistical accumulation 
of all possible outcomes in order to faithfully replace frequency of outcomes with 




Figure 5. (a) Plot showing the twofold probabilities using a |1,1) input state 
propagated through three levels, (b) Feynman diagram showing one set of unique 
paths that reach detectors Dl and D5. (c) Feynman diagram showing the only other 
unique set of paths that the photons can take to reach detectors Dl and D5. 



probabilities. In particular, detection at I level has 21 outputs. Therefore for the N- 
walker case there are (2/)^ outcomes that carry information about the quantum random 
walk. Each particular outcome provides only limited information about the quantum 
walk. In order to get an intuitive understanding of what is happening to the photons 
in quantum random walks, we use Feynman diagrams. We found that it is important 
to grasp how the phase shifts are being picked up and how they create constructive and 
destructive interference, visually. In the following sections, we will show specifically how 
the diagrams work, using examples of input states with a given number of photons. 



4-2. Onefold probability 



The complete destructive interference at detector D3 shown in Figure 4(a) distinguishes 



it from the classical random walk and can be explained using a Feynman diagram. 



Figure 4(b) and Figure 4(c) show the Feynman diagrams for the onefold probability of 
detector D3 firing with the |1, 0) input state after three steps. These diagrams show all 
possible paths for the single photon to reach detector D3. In this case there are only 
two possible paths, which means that no matter what is happening to the photon on all 
other paths, this particular outcome is unaffected. 

In order to calculate the destructive interference, we define t and r as the coefficients 
that a photon's probability amplitude accumulates when it is transmitted or reflected. 

t = l/V2, (Photon is transmitted), 

r = i/v / 2, (Photon is reflected). 

The definitions of t and r remain consistent with the beam-splitter transformation ([!} 



(3) 



When the photon follows path A in Figure 4(b) , it is reflected off the first, second, 
and third beam-splitters. Therefore this path has a total probability amplitude of 
PathA = r 3 = — i/2 3//2 . When the photon follows path B in Figure 4(c), it is 
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(c) 

Figure 6. Plot of threefold probabilities using input state of |2, 1) propagated through 
3 levels of beam-splitters. Each row, column, and pillar represents an individual 
detector while the color of each cube represents the probability. Blue represents 
minimum probability, while red represents maximum, and white is exactly zero, (a, 
b) Plots of the threefold probabilities using input state of |2,1) propagated through 
3 levels of beam-splitters detection scheme looking into the space diagonals. These 
views display the symmetry along the diagonals, (c) "Exploded" view of the same 
parameters. This view allows you to see inside the cube. 



transmitted through the first beam-splitter, is reflected off the second, and is transmitted 
through the third beam-splitter. Therefore as the photon follows path B, it receives 
a total probability amplitude of PathB = trt = i/2 3//2 . In order to calculate the 
probability of photon detection at D3, the probability amplitudes are added together, 
PathA + PathB = -i/2 3/2 + i/2 3/2 = 0. All the paths combine to zero, which confirms 
our result in Figure 4(a) at detector D3. 

Constructive interference at D4 is revealed by looking at the Feynman diagrams in 
Figure 4(d) and Figure 4(e) , when we follow path A, we see that there are two reflections 
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Figure 7. Feynman diagrams explaining the threefold probability at detectors D3, 
D4, and D6 using a |2, 1) input state. Each diagram is a set of unique paths for 
detectors Dl, D4, and D6 to fire simultaneously along with their respective probability 
amplitudes. We have each diagram as a separate possible case of photon paths. 
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followed by a transmission. This gives us PathA = r 2 t = — 1/2 3 / 2 . Following path B, we 
see one transmission followed by two reflections. This gives us PathB = tr 2 = — 1/2 3 / 2 . 
In order to calculate the probability of photon detection at D4, we add these two 
probability amplitudes PathA + PathB = — 1/ y/2. This is the probability amplitude of 
the photon being detected at detector D4. So then if we look at the probability | — 1/ \/2\ 2 
we see that the probability of the photon reaching detector D4 confirmed by our 



simulation in Figure 4(a) 



D3 and D4 present information about essentially the same path in the interferometer 
except for the final leg. However, from an experimental perspective, D3 is more 
informative since seeing a photon there tells us that something happened to the system. 
However, detecting or not detecting a photon at D4 does not immediately tell us that 
something has changed until one accumulates enough statistics. 



4-3. Twofold probability 



In order to obtain more information about the system, one needs to go to the higher order 
correlation function and multi walker cases. Let us consider destructive interference 
with two walkers. Using the t and r operators defined above fl3]) and referring to the 



Feynman diagrams in Figure 5(b) and Figure 5(c), we calculate the twofold probability 



of detectors Dl and D5 firing simultaneously. Figure 5(b) shows one possibility of the 



two photons reaching detectors Dl and D5. These paths are described by 



Figure 5(b) 



PathA 
PathB 



[rt 2 ) 
(rtr) 



i/2 3 / 2 
= -1/2 3 / 2 



(4) 



Since these paths take place at the same time, we multiply these amplitudes 



together, (i/2 3//2 )(— 1/2 3 / 2 ) = —i/2 3 . Figure 5(c) shows the only other paths the photons 
could take in order for detectors Dl and D5 to fire. These are described by 



Figure 5(c 



PathC 
PathD 



{t 2 r) 
(t 3 )- 



= i/2 3 / 2 

1/2 3/2 



(5) 



Similarly, these two paths cannot take place independently of each other, so we 
multiply them, (i/2 3 / 2 )(l/2 3 / 2 ) = i/2 3 . Adding two possible scenarios is necessary in 
order to calculate the probability of a joint event, (i/2 3 / 2 ) (— 1/2 3 / 2 ) + (i/2 3 / 2 )(l/2 3 / 2 ) = 
— i/2 3 + i/2 3 = 0. This complete destructive interference is shown in our plot of 
the twofold probability at detectors Dl and D5 in Figure 5(a) One cannot predict 
this destructive interference from the onefold probability. It truly is multiphoton 
interference! 



4-4- Threefold probability 

We again use Feynman diagrams to help us explain the values of the threefold 
probabilities. This detection scheme is more complicated than the previous detection 
schemes and thus the Feynman diagrams are also more complex. These diagrams become 
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more complex because we have at least three photons in this setup but only two possible 
input ports. This introduces normalization factors that were not present in our previous 



examples. Looking at Figure 7(a) we see one set of unique paths for threefold coincidence 
detection with an input state of |2, 1) and their associated probability amplitudes. 
Analyzing this diagram further we have 

i/2 3 / 2 
-1/2 3 / 2 . 



frr 



Figure 7 (a 



PathA 
PathB 
PathC = 1 2 



(r 2 t) 



(6) 



i/2 3 / 2 



Since each of these paths are independent, we multiply their individual probability 
amplitudes. 

PathA ■ PathB ■ PathD = (i/2 3/2 )(-l/2 3/2 )(i/2 3/2 ) = l/2 9/2 . 



Here we note that Figure 7(b) is the same as Figure 7(a) except the photons that share 



an input port are switched. Since this switching does not affect the output distributions 
we will get the same probability amplitude of l/2 9//2 for this case. 

We perform the same analysis for the remaining Feynman diagrams. For each 
diagram we simply take the product of the probability amplitudes as we did before. 

1/2 9 / 2 
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(7) 



Now to get the total probability of detectors Dl, D4, and D6 firing simultaneously 
we take the sum of each of these twelve cases. 

Total Probability Amplitude = 6(l/2 9 / 2 ) - 6(l/2 9 / 2 ) = (8) 
So we expect to see complete destructive interference for the threefold probability 



involving detectors Dl, D4, and D6. This is shown in Figure 6(c) as block (1,4,6) 
(4, 1, 6) = (6, 1, 4) = (1, 6, 4) = (4, 6, 1) = (6, 4, 1) 



4-5. Rules for construction of Feynman diagrams 

Feynman diagrams allow us to describe specific scenarios within our QRW in terms of 
intuitive quantum interference. It is advantageous then to have a set of rules dictating 



QRW with multiphoton interference and high order correlation functions 



13 



how to generally construct these diagrams, based on the choice of, input state, detector 
combination, and level of propagation. 

The number of paths in each Feynman diagram is equal to the total number of 
photons in the system. The use of multiple diagrams for a single choice of input is then 
due to the fact that each photon is indistinguishable from another; therefore we must 
account for every photon reaching each of the desired detectors as separate cases. Using 
Figure 5(b) and Figure 5(c) as examples, there are two paths per diagram and two 
diagrams because there are two photons in this scheme but either photon can impact 
either of the detectors. 

We define a pair of detectors as any two detectors that share the output of a common 
beam-splitter. The number of paths that lead to any pair of detectors are given by the 
binomial expansion (x + y)^ 1-1 ' where I is the level of propagation and the coefficient 
of the i th term gives the number of paths for the i th pair of detectors. This shows that 
pairs of detectors that reside on the edge of the system only ever have one path that 
lead to them, regardless of level. It also shows that any non-edge detectors, which we 



call central-pair detectors, have multiple paths that lead to them. Using Figure 7(a) 



and Figure 7(c) as an example, detectors D3 and D4 form a central pair and each have 
two paths that lead to them. 

Another case that leads to multiple diagrams is when multiple photons are 
introduced from the same input port. This is always the case in the threefold detection 
scheme, as there are always three total photons but only two input ports. As before, 
since all the photons are indistinguishable from one another, we must have a diagram 
for each photon's path to each detector, regardless of if it emerges from the same input 



port as another photon. An example of this is shown in Figure 7(a) and Figure 7(b) 



Notice that swapping the paths of the two photons (red and blue path) that share an 
input, does not change the output distribution. 

With these rules, in principle, one could use Feynman diagrams to describe any 
scenario within our system, but realistically, we limit the use of Feynman diagrams to 
describe low-level scenarios to reduce the number of required diagrams. The number of 
Feynman diagrams grows exponentially with increasing photon number and level as the 
associated Hilbert space needed to describe the system also grows exponentially with 
these parameters. 



5. Future research 



5.1. Twofold number resolving detection 

With the recent development and use of number-resolving photon detectors, we propose 
a new detection scheme that can be utilized [14]. Our twofold-number resolving 
detection scheme is an extension of the twofold-detection scheme but with the added 
benefit of number resolving detectors. In this detection scheme we can investigate the 
probability of detector m receiving exactly i photons given detector n receives exactly j 
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photons simultaneously. The result of this detection scheme would be four dimensional 
and would be much more difficult to visualize. We are currently analyzing the results 
of this detection scheme. The benefit of this scheme over the standard twofold scheme 
are the cases where there are more than two photons present in the system. 

We are interested in investigating further modifications to the cascading beam- 
splitter structure. In our current simulation, we only introduce vacuum into the edge 
ports at each level. Introducing photons into the side ports in a systematic way may 
result in more interesting and useful output states. Another avenue could be the use 
of phase shifters in between specific beam-splitters. We could change the amount of 
relative phase shift as well as the placement of these phase shifters. If placed properly, 
these phase shifters may effectively steer photons to a desired port with constructive 
interference or create a dark port by canceling out all phases at that detector. The 
possible combinations of placement and phase shifts would surely be enormous so we 
have considered implementing a genetic algorithm or similar technique to fine tune the 
system in order to create a desired output. 

6. Conclusion 

We have shown that quantum random walks exhibit properties that are not present 
in their classical counterpart. Furthermore, multiphoton interference gives rise to new 
phenomena not found in standard quantum random walks. These properties are due to 
the quantum nature of the input photons and gives rise to unique interference effects. We 
simulate this quantum random walk with an array of 50/50 beam-splitters and a bank of 
detectors at any desired level of the system, only limited by the computational resources 
required to calculate these schemes. It has been suggested that these random walks have 
applications in quantum computing and solving the graph isomorphism problem [El [THE]. 

We expanded upon the onefold probability detection scheme with our two and three 
fold coincidence detection schemes. Two fold detection has been studied previously 
but our three fold detection scheme, while harder to visualize, is the next logical step 
in simulating a quantum random walk and provides more information than a twofold 
probability scheme [1]. 

We use Feynman diagrams to give an intuitive explanation of how the interference 
effects propagate through our system and result in the final output distribution. These 
diagrams are functionally limited to low levels of propagation but still act as a valuable 
tool in order to explain how our simulation propagates the input photons. 
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